Draft VERSION April 25, 2013 

Preprint typeset using MrgX style emulateapj v. 5/2/1 1 



STUDYING INTER-CLUSTER GALAXY FILAMENTS THROUGH STACKING GMBCG GALAXY CLUSTER PAIRS 

YuANYUAN Zhang', JorgP. Dietrich'"\ Timothy A. McKay'-^-'*, and Erin S. SHELDON^ AlexT. Q. Nguyen''' 

Draft version April 25, 2013 



m 

o 

(N 

Vh 

< 

(N 

o 
u 

p.; 
o 



> 

o 



% 



ABSTRACT 

We present a method to study the photometric properties of galaxies in filaments by stacking the galaxy 
populations between pairs of galaxy clusters. Using Sloan Digital Sky Survey data, this method can detect the 
inter-cluster filament galaxy overdensity with a significance of ~ 5a out to z = 0.40. Using this approach, we 
study the g-r color and luminosity distribution of filament galaxies as a function of redshift. Consistent with 
expectation, filament galaxies are bimodal in their color distribution and contain a larger blue galaxy population 
than clusters. Filament galaxies are also generally fainter than cluster galaxies. More interestingly, the observed 
filament population shows redshift evolution at 0.12 < z < 0.40; the blue galaxy fraction increases at higher 
redshift: a filament "Butcher Oemler Effect". We test the dependence of the observed filament density on the 
richness of the cluster pair: richer clusters are connected by higher density filaments. We also test the spatial 
dependence of filament galaxy overdensity: this quantity decreases when moving away from the inter-cluster 
axis between a cluster pair This method provides an economical way to probe the photometric properties of 
filament galaxies and should prove useful for upcoming projects like the Dark Energy Survey. 
Subject headings: Large Scale Structures: Galaxy Filaments: Photometric Properties 



1. INTRODUCTION 

In the standard picture of structure formation after the 
Big Bang, initial density perturbations grow via gravita- 
tional instability and form massive structures in a bottom- 
up, hierarchical growth scenario in which small, gravita- 
tionally bound structure like galaxies conglomerate into pro- 
gressively larger structures like galaxy groups and clusters. 
On the largest scales, gravitation al instability lead s to a col- 
lapse of matter first into sh eets (|Zeldovichlll970l) and then 
into a network of filaments (lKlvpin& Shandarinlll983h with 
galaxy c lusters at the intersec tion of filaments. This early 
result of iKlypin & ShandannI has been c onfirmed in count- 
less A^-body simulations since then (e.g.. | Davis et alJ 119851 : 
iBertschinger & Gelbl[T99Tt ISpringel et al.ll2005l) . 

Observationally the filamentary large-scale structure (LSS) 
has been traced by the galaxy distribution in redshift sur- 
veys f or decades; from first indic ations seen by H Beveer et al] 
(Il978h over the seminal work of iGeller & Huchra (1989) to 
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SDSS ( lAhnetal.1120121) . The gaseous Warm-Hot Intergalac- 
tic Medium residing in filaments was se en in X-ray emis- 
sion (IWerneretal.ll2008l) and absorption (iBuote et al.l 120091 
iFang et al.ll2010h . More recently the underlying Dark Matter 
skeleton of large-scale structure filame nts was also detected 
through its gravita tional lensing effect (iDietrich et al.l 120121 : 
iJauzac et alj|2012l) . 

Despite these observational advances, filaments remain rel- 
atively little studied, yet they are astrophysically interesting 
places: LSS filaments contain a plura lity of all matter in 
the Universe (lAragon-Calvo et aDl2010l) and they harbor the 
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"missing baryons" at low redshift (iDave et al.l2001l) . It is well 
known that galaxies in clusters are redder than galaxies in the 
surrounding field and mostly have ceased star-formation. As 
galaxies are constantly accreted into clusters along filaments, 
it is only reasonable to assume that a portion of the transfor- 
mation from actively star-forming galaxies to passively evolv- 
ing ones happens in filaments. The specific role of the fil- 
amentary environment in this process has not been studied 
very m uch. An exception is, e.g., the work of iPorter et al.l 
(12008). which found a spike in star-formation activity in fila- 
ment galaxies at fixed cluster centric distance. It is of course 
well established that the local density of galaxies is one of 
the factors d eciding the efficie ncy with which star-formation 
is quenched (iPenget alJl2010l) . The exact role of the filament 
versus the galaxy-group environment, however, has not been 
studied in detail. 

The relative dearth of papers about filaments compared to, 
e.g., the number of studies done of galaxy clusters is due to 
the low-density contrast of filaments with respect to the mean 
density of the Universe. While galaxy clusters by definition 
have a density contra st 5 > 200, the typical den sity contrast 
of filaments is (5 < 20 dAragon-Calvo et al.l2010l) . This makes 
filaments much more difficult to find, observe, and study, and 
they have generally been traceable in spectroscopic surveys 
only to redshift ~ 0.2. 

In this paper, we develop an algorithm to study the photo- 
metric properties of inter-cluster galaxy filaments from red- 
shift 0.1 to 0.4. Because close pai rs of galaxy clusters 
are generally conne cted by filaments (iPimbblet et al.l l2004t 
IColberg et al.l l2005h . instead of trying to detect single fila- 
ments, our approach is to identify galaxy cluster pairs and 
then stack the galaxy population between them. By compar- 
ing the stack of such populations to the stack of other selected 
galaxy fields, we can detect an overdensity of galaxies present 
in the connecting filaments. We go on to analyze the color and 
luminosity properties of these filament galaxies and to exam- 
ine these for possible redshift evolution. 

The remainder of this paper is arranged as follows. In Sec- 
tion 111 we describe our algorithm and some tests we have 



conducted on it, and represent the statistical significance of 
our filament overdensity. Section |3]contains our results. We 
compute color evolution and luminosity function of filament 
galaxies, as well as their spatial and richness dependence. We 
summarize and discuss the implications of our results in Sec- 
tionlH In AppendixlAl similar results from N-body simulation 
are provided for comparison. Throughout this paper, we as- 
sume a flat cosmology with fJn, = 0.30 and h = 0.7. We use z 
to denote redshift, and S to denote the 2-d projected galaxy 
number density. All the error bars shown in this paper repre- 
sent single standard deviation errors. 

2. METHODS 

2.1. Data 

For th is study, we use the gmBCG galaxy cluster catalog 
(iHao et al. 2010). This catalog is const ructed from the SDSS 
Data Release 7 (lAbazajian et al.ll2009l) . and includes 55,424 
clusters in the r edshift ran ge 0.1 < z < 0.55. When build- 
ing this catalog, iHao et all searched for overdensities of red 
sequence galaxies around brightest cluster galaxy (BCG) can- 
didates. The final catalog contains positions, estimated red- 
shifts, and richness for each detected cluster 

The galaxy catalog we use for thi s study is constructed 
from SDSS DRS B OSS imaging data (lEisenstein et alJl2011l: 
lAihara et al.ll201 ll) . The sky coverage of this catalog is de- 
scribed by a catalog of random points, which samples the 
same angular coverage as the survey objects. Instead of us- 
ing SDSS DR7, we choose to use galaxy data from the BOSS 
survey in a newer SDSS data release because of its improved 
photometry. This choice leaves the galaxy catalog having a 
slightly different sky coverage with the gmBCG cluster cata- 
log. H owever , with the sky coverage weighting described in 
Section l2.2.5l only the overlapping area of the galaxy catalog 
and the cluster catalog contributes to our final results. 

2.2. Algorithm 

Figure [T] outlines our algorithm, which is described in de- 
tail in this section. To summarize, we select pairs of galaxy 
clusters and identify the galaxies which occupy the regions 
around these cluster pairs. We stack these galaxies, weight 
for sky coverage, and compare the galaxy population between 
cluster pairs to the galaxy population seen in other fields. 

To test our methods and compare to our SDSS results, 
we also apply our a lgorithm to the Millennium simulation 
(ISpringel et al.ll2005h . The procedures and results of this ap- 
plication are discussed in the appendix. 

2.2.1. Cluster Pair Selection 

We begin by selectin g the cluster pairs to b e used in our al- 
gorithm. According to IColberg et al.l (120051) . the probability 
of finding a filament between galaxy clusters is strongly de- 
pendent on their separation. Clusters separated by less than 
5/!~'Mpc always have a connection between them, but this is 
primarily formed from the outer regions of the cluster popu- 
lations rather than a separate filament population. For clus- 
ters that are separated between 15/i~'Mpc to 25/!~'Mpc, about 
1 /3 of them have a filament connection. For clusters that are 
separated less than 35/i~'Mpc, the probability of filament ex- 
istence between them is larger than 10%. Clusters separated 
more than 50/!~'Mpc are unlikely to have filaments between 
them. Also, richer clusters tend to have more filaments con- 
necting to them, but the number of connections rarely exceeds 
5. Informed by these expectations from simulation, we select 



cluster pairs separated by less than ^ 35/r'Mpc while lim- 
iting the maximum number of cluster pairs one cluster can 
appear in. For each cluster, the clusters paired to it satisfy the 
following conditions, 

1 . The difference in photometric redshift of the BCGs of 
the two clusters should differ by < 0.02. This redshift 
difference roughly corresponds to a comoving distance 
of50/riMpc. 

2. The angular separation between the two clusters is 
smaller than 35/r'Mpc on the plane of the sky at the 
redshift of the the cluster to-be-paired-to. Because we 
do not directly constrain the two clusters' 3-d physi- 
cal separation, the cluster pairs' 3-d physical separation 
can be larger than 35/i~'Mpc, as shown in Figure|2] 

3. Each cluster can appear in several cluster pairs if it has 
more than one neighbor meeting the above criteria. If a 
cluster appears in more than 5 cluster pairs, we keep its 
5 closest pairs along with all pairs that include clusters 
with fewer than 5 neighbors. 

With this selection procedure, we ide ntify a total of 160,954 
cluster pairs. As discussed in Section 12.2.61 when stacking 
galaxy populations that are around these cluster pairs, we ap- 
ply an absolute magnitude cut which is above the SDSS com- 
pleteness limit at z = 0.38, even with dust extinction correction 
and K-correction. To avoid including galaxy populations in 
very dusty sky regions, we exclude cluster pairs which fall 
within 1° X 1° angular boxes around galaxy objects which 
have dust extinction in r-band larger than 0.4. This lowers 
the number of cluster pairs to 158,897 in the final catalog. In 
this cluster pair catalog, one cluster on average appears in 5.7 
cluster pairs. 

2.2.2. Selection of Galaxies 

For every cluster pair, we identify all galaxies in and around 
the two clusters. We first identify all galaxies in a rectangular 
box centered on the filament midpoint, with the width of the 
inter-cluster distance but twice as long. Galaxies in the square 
centered on each cluster form the "galaxy set" for that cluster 
Each galaxy set includes the cluster, an inter-cluster region 
in which we expect to find a filament, and an opposing field 
region in which we do not (necessarily) expect a filament. 

When selecting these neighboring galaxies, we do not make 
any redshift cuts. It is possible that a carefully designed pho- 
tometric redshift cut coul d inc rease the significance of the re- 
sults discussed in Section l231 but such a cut risks to introduc- 
ing color and magnitude dependent selection effects. While 
photometric redshifts fo r bright passively e volving galaxies 
can be smaller than 0.05 (ICunha et al.ll2009h . those for fainter 
or star-forming galaxies can be substantially worse. 

2.2.3. Cluster Galaxy Set Processing 

For each cluster in one cluster pair, there is one galaxy set 
corresponding to it. We denote the angular coordinates of one 
cluster in one cluster pair to be (ao, So), its comoving dis- 
tance to be dc, the angular coordinates of the ^th galaxy in 
its galaxy set to be {a„,6s), and the angular coordinates of the 
other cluster BCG in this cluster pair to be (aop,(5op)- Before 
stacking the galaxy sets of different clusters, we need to rede- 
fine galaxy coordinates so that they convey the same physical 
scale at different redshifts; 
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Fig. 1. — Flowchart illustration of the algorithm described in Section [Z2l The left side of the figure illustrates the manipulations on galaxy data in order 
to acquire stacked cluster pair galaxy fields. The right side illustrates the the manipulations on the random point catalog in order to acquire the sky coverage 
weighting, which tells how many times each pixel in the galaxy stacking is covered by real sky data. The manipulations on random points are identical to the 
manipulations on real galaxy data before "Compute the Weighting" step (See Section 12. 2. SI for details). Aitows indicate the processing flow. We first look for 
cluster pairs that possibly have filament connections and select the galaxies that are in the fields of such pairs. We then re-define the coordinates of thes e galax ies 
with steps "Project", "Translate" and "Rotate" so that the coordinates reflect the galaxies' projected physical distance to cluster centers ( See Section |2. 2. 3| for 
details). We stack these galaxies and weight the stack with sky coverage weighting to obtain "averaged" galaxy distributions around one cluster center in cluster 
pairs. See the electronic edition of the journal for a color version of this figure. 
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Fig. 2. — Cluster pairs' projected separation on the plane of the sky at their 
median redshift, rp,oje;.tcd. versus their 3-d physical separation, r^o- The white 
dashed line shows the mean rjQ value at different ^projected- We assume the 
clusters' photometric redshift are purely due to Hubble flow and use the law 
of cosine for calculating r^o- The cluster pairs' projected separation on the 
plane of the sky is cut off at ~ 35/r ' Mpc, depending on the exact redshift of 
the projection, but their 3-d separation can be as large as 90/;"' Mpc. 

1 . Project these galaxies onto a tangent plane at the cluster 
BCG, 



9g = agCos{So), 



(1) 



so that 1 ° in the 6^ and Sg dimension corresponds to the 
same great-circle distance. 



2. Translate {Oa,Sg) as in 

e'g = eg-(aocos(So)), 
5g = 5g-So, 



(2) 



so that the cluster center of the galaxy set is placed at 
(0,0). 



3. Rotate the galaxy set with a rotation matrix M, 



Si' 



■ Mx 



(3) 



to have the inter-cluster axis lay on y = 0, with the direc- 
tion extending to the other cluster aligned toward x > 0. 
The rotation matrix M has the form of: 
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4. Redefine the coordinates to be 

x=6'^'x(c/ctanl°), 
y = S'J x{dctiinl°). 



(4) 



(5) 



Here, dc tan 1 ° is the physical distance that corresponds 
to 1 ° angular separation on the tangent plane at the clus- 
ter center If all the galaxies in the galaxy set are at 
about the same redshift with the cluster BCG, the co- 
ordinates re-defined above would reflect the projected 
physical distance between these galaxies to the cluster 
BCG in the plane of the sky. This argument does not 



hold for the whole population since we do not make any 
redshift selection. However, this definition of coordi- 
nates should be proper for filament and cluster galaxies 
after foreground and background subtraction. 

2.2.4. Stacking 

When stacking, we create a multidimensional galaxy stack 
in bins of cluster pair redshift. Each stack is four dimen- 
sional, including galaxy coordinates {x,y), galaxy absolute 
magnitude M,-, and the observer frame dust extinction cor- 
rected g- r color The binning size in each dimension is 
chosen to be small enough not to smear out the details of 
X, y, Mr and g-r dependence, but not so small as to leave 
many bins unoccupied. For the x and y dimensions, the bin 
sizes are both 0.7/r'Mpc, because of the low S/N in one 
spatial bin. For binning in Mr and g-r, the bin sizes are 
0.2 mag and 0.125 mag. The value of the stack at one partic- 
ular bin Xi-i <x< Xi, Yj-\ <y < Yj, Mru-\) <Mr < Mri and 
(8-r)k-i < (g-r) < {g-r)k is denoted as G;(x/,yy,Mw,(g- 
r)k). To include a comparison of cluster galaxies to fila- 
ment galaxies, we also build a two dimensional stack in Mr 
and g-r with galaxies around cluster centers which satisfy 
yjx^+y^ < 0.7/;"' Mpc. This stack of cluster galaxies is de- 
noted as CL^ (Mr ,g-r). 

2.2.5. Stack Weighting 

Because the galaxy cluster pairs have different separations, 
and the galaxy sets of different clusters are selected based on 
their cluster pair separations, the physical sizes of our clus- 
ter galaxy sets vary. Also, the cluster pairs may fall at the 
edge of the survey, and the galaxy sets will have irregular sky 
coverage. Without accounting for this, we would observe ar- 
tificial galaxy density gradients in the x and y dimension of 
the stacking. To remove such an artificial gradient, we weight 
G^{x,y,Mr,g-r) with a coverage weighting function W7{x,y). 
We first explain the physical interpretation of W~{x,y) and the 
weighting procedure in this and the next paragraph, and de- 
scribe details on computing Wr{x,y) in the rest of this section. 
The value of W~(x,y) at X,_i < x < X, and Yj-i <y< Yj, de- 
noted as W-{Xi,Yj), tells how many cluster galaxy sets include 
valid data for the pixel at X,_i < x < X, and Yj-i < y < Yj. 
For example, W~(Xj,Yj)=3 when there are 3 clusters' galaxy 
sets fully covering this pixel, or W-(Xi,Yj) = y when there are 
3 clusters' galaxy sets fully covering this pixel and another 
cluster's galaxy set covering 1/3 of this pixel. 

We weight G-(x,y,Mr,g-r) with Equation|6] and acquire a 
weighted distribution, g-{x,y,Mr,g-r), the value of which at 
one bin, Xi-i <x < Xi, Yj-i <y < Yj, Mr(i-i) < Mr < Mri and 
(8-r)k-i <i8-r)< ig-r)^, is denoted as g-Xxi,yj,Mri,{g- 
r)&. 



8z(xi,yj,Mri,(8-i-)k) = 



G:(Xi,yj,Mrl,(8-l')k) 

W,(xi,yj) 



(6) 



This weighted distribution, g^(x,y,Mr,g- r), can be inter- 
preted as representing properties of one "averaged" cluster's 
galaxy set with filament connection on the x > side. 

The W~{x,y) we u se is built from the random point catalog 
described in Section lTTI This random point catalog first goes 
through the same algorithm with t he gal axy catalog as de- 
scribed in Section l2.2.2l a nd Sec tion l?.2.3l and then is stacked 
together Unlike Section 12.2.41 in which we build a four di- 
mensional distribution from stacking galaxies, we can only 



build a two dimensional distribution in x and y from stack- 
ing random points, denoted as R,{x,y). At X,_i < x < X, and 
Yj-i <y< Yj, the value oiR-Sx.y), denoted as R-XXt^Yj), tells 
how many random points there are in this bin. Because the 
random point catalog covers the survey area evenly, /?-.(x,y) 
can be used as a proxy of W.<(X,y). The relation between 
R-{x,y) and WJ^x^y) at one particular bin is 



R,{xi,yj) = W,(,Xi,yj)xC{z). 



(7) 



Here C(z) is the average number of random points per 
[/i"'Mpc]^ at redshift z, and z is the median redshift value of 
the clusters which have their galaxy sets or random point sets 
stacked to build G~{x,y,M-f) or Rr{x,y). C{z) can be com- 
puted from the comoving distance, denoted by dciz), at red- 
shift z, and the average number of random points per deg , 
denoted by C, through 



C{z) = C X 



arctan 



0.7/i-'Mpc 
ddz) 



(8) 



With known values of RXx,y), and C, We compute WXx,y) 
from Equation |2] and Equation |8] and then gJ^x^y^Mi^g- r) 
from Equation |6] which is used for deriving the scientific re- 
sults in this paper. 

To weight CL.{Mr,g- r), which is used for comparison 
with g-{x,y,Mr,g- r), we record the total number of ran- 
dom points that fall around cluster centers with i/x^+y^ < 
0.7/2"'Mpc, denoted by R-, and compute the weighted func- 
tion of CL~{Mr,g-r), which is denoted by cL(Mr,g-r), with 
following equations: 



ch{Mr,,(g-r)k) = 



R, = W,x C{z), 

CL,(Mri,(8-r)k) 



(9) 



W^ 



In this weighting procedure, one might worry that the pois- 
son noise of the random point sampling would lower the S/N 
of g{x,y, Mr, g-r) or cl.{Mr,g- r). Here, we show that with 
proper random point density, this influence is insignificant. 
The poisson noise in R-{x,y) and R~ increases with redshifts 
because C{z) decreases with redshifts. Our random point cat- 
alog samples the survey coverage with C = 9,275 and at z = 
0.40, C{z) = 5.98. Since we always stack galaxy sets of more 
than 5,000 clusters in each redshift bin, the S/N of R-{x,y) at 
most spatial bins, and of/?-, is larger than 100. At this signif- 
icance level, the weighting constructed from R~{x,y) and R~ 
brings in negligible noise to g{x,y,Mr,g-r) or cl,(Mr,g-r). 

2.2.6. Galaxy Absolute Magnitude and Luminosity Cut 

For every galaxy selected in Section 12.2.21 we compute an 
absolute magnitude from its apparent magnitude, nir, SDSS 
dust extinction correction, e, K-correction, K, and its cluster 
BCG's luminosity distance d]^, with equation 



Mr = nir - 5 log 



lOpc 



-e-K . 



(10) 



We compute the K-correction fo r each galaxy with analyt- 
ical approximations provided by iChiUnparian et al.l (120101) . 
since the more popular SED template fitting method 
is inefficie nt to implement with our large volume of 
data. This IChiUngarian et al.| approach ap proximates K- 
corrections computed by iFioc & Rocca-Volmerange, (,1997,) 



and lBlanton & RoweisI ( 120071) with analytical polynomials of 
two parameters: redshift and one observed color, and the 
residual between values measured by this approach and others 
fepc & R occa-Volmerangd 119971: [glanton & RoweisI 120071: 
iRoche et al.112009 ') is clos e to zero for r-b and in the redshift 
range [0.1,0.45] (Chiling arian et al.ll201CII) . In our application 
of this method, we use each galaxies own g-r color and the 
cluster BCG's photometric redshift. 

Because SPSS is an apparent magnitude limited survey 
dStoughton et al.l 120021) . its imaging data contains different 
galaxy populations at different redshifts: at lower redshift, 
the galaxy catalog is more complete, including more intrin- 
sically faint galaxies than at higher redshift. To ensure we are 
selecting the same galaxy population for comparison in differ- 
ent redshift bins, we apply a -24 .2 < M, <-20.4 cut in Sec- 
tion|231 Section|23] SectionO Section[33]and Section[14l 
Under this magnitude cut, even with a maximum dust ex- 
tinction correction, Cmax = 0.4, and a maximum K-correction 
^max = 0.6 in the galaxy catalog, a galaxy's absolute magni- 
tude is still above the 95% completeness limit of r-band at 
redshift z = 0.38. In Section [32] the M,. cut is different for 
different redshift bins, which is listed in Table [1] This cut en- 
sures that the galaxy population is more than 95% complete 
in each redshift bin. 

2.3. Foreground/Background Subtraction and Filament 
Detection Significance 

In this section, we describe our foreground and background 
galaxy subtraction and quantify the detection significance of 
our filament signal. Within the weighted average cluster 
galaxy set g.{x,y,Mr,g-r) we define a "filament" region (de- 
noted FL), a "cluster" region (denoted CL), and two compari- 
son regions CI and C2. All of these regions are displayed in 
Figure [3] We search for filament galaxies by subtracting the 
galaxy population observed in CI or C2 from that seen in FL. 

Since the filament fields are placed on the right side of the 
galaxy sets, the "filament" region is defined with the following 
constraints: 

2. 1 /r' Mpc < X < 10.5 /j"' Mpc 
|6» = arctan(3;/x)|<45°. 

We define this region to be 2.1/1"' Mpc away from the cluster 
center to avoid including contamination by cluster galaxies. 
Although cluster galaxy contents can be removed with the 
foreground/background subtraction procedure as described 
later in this section, contamination from this population would 
lower the S/N of filament signals. We also eliminate galax- 
ies that are too far away from cluster centers (with x > 
10.5/2~'Mpc), to avoid lowering the S/N of filament detec- 
tions with too noisy spatial bins. Since few cluster pairs 
extend to 10.5/r'Mpc, at the spatial bins beyond this dis- 
tance, not only the galaxy count noise becomes significant, 
the random point count noise also starts to influence, result- 
ing in exceptionally low S/N of these bins. The filament re- 
gion angu lar e xtent 16*1 = 45° is chosen based on discussion 
in Section 13.41 Defining the cluster region is more straight- 
forward. Since cluster BCGs are placed at {x,y) = (0,0), we 

define the region enclosed within y^x-+y^ < 0.7/!"'Mpc to 
be the cluster region, which is also the coverage of functions 
CL,{Mr,g-r) and cL{Mr,g-r). 

To evaluate signals from the foreground/background galaxy 
population, we define two comparison regions. The first, CI, 
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Fig. 3. — (a) Shows a galaxy overdensity significance contour map from stacking galaxy sets of clusters at 0.14 <z< 0.18. There is a significant overdensity 
in the middle of the map, caused by the presence of many cluster galaxies, (b) Shows null test results from stacking randomly re-positioned pairs, where the 
whole field is noisily flat. In (a) and (b), the solid fine and dashed line boxes mark out the four regions defined in Section 1231 See the electronic edition of the 
journal for a color version of this figure. 



satisfies these constraints: 

-10.5/i"'Mpc <x< -2.1/r'Mpc 
|6' = arctan(3;/jc)| <45°, 

The second comparison region, C2, satisfies, 

-13.3/r'Mpc < X < -6.3/r'Mpc 
-10.5/i"'Mpc <3/ < 10.5/2"'Mpc. 

The CI region is symmetrical to the FL region, and is used 
for filament foreground/background subtraction. Both FL and 
CI contain galaxies present in the outskirts of the clusters, 
so using CI for foreground/background subtraction should, to 
first order, eliminate these. The C2 region is used for fore- 
ground/background signal evaluation in estimating the cluster 
population. It is picked to be a region away from the cluster 
center 

In Figure m we resample with replacement, i.e. bootstrap, 
cluster pairs which have their galaxy sets stacked, and show 
the projected galaxy number density in the FL region and the 
CI region. The averaged galaxy number density in the fila- 
ment region is ~ 5a higher than that of the CI region. Since 
cluster outskirt galaxies have already been taken into consid- 
eration when making comparison, this galaxy overdensity is 
unlikely to be caused by the existence of one galaxy cluster 
alone, but related to cluster pair structures exclusively. We 
associate this galaxy overdensity with the filament population 
between cluster pairs. Such a galaxy overdensity is very small 
comparing to galaxy overdensity caused by cluster fields. In 



Figure lUa), we show projected galaxy overdensity signifi- 
cance contour map of stacked galaxy sets at 0.14 < z < 0.18. 
At {Xi^Yj), the projected galaxy overdensity significance, de- 
noted by S{Xi,Yj), is defined as. 



S,{X;,Yj) = 



Y.gzi^i^yiAM,)tAg-r)i)-tiz 

H (11) 



Here, fj,- and ct,, are the mean and scatter of projected galaxy 
count per spatial pixel in the C2 region, computed with. 



fJ-z = 



Y: 'E8z(Xi,YjAM,-)kA8-r),) 

a J) in C2 kj 

(ij) in C2 



a- = 



\ 



(iJ) in C2 k.l 

(iJ) in C2 



(12) 



In Figure |3ja). We observe a significant overdensity in the 
center because of cluster galaxy population. The filament 
overdensity is not directly observable in this figure, due to 
the very small density contrast of the filaments. 

Because more than 90% of the galaxies observed in the FL 
region come from foreground and background populations, 
the properties of filament galaxies would not be distinguish- 
able without foreground/background subtraction. To get the 
pure properties of galaxy filaments, like color distribution and 
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Fig. 4. — Projected galaxy number density in the filament region (FL) an d the comparison 1 region (CI) at 0.14 < z < 0.18 when stacking real cluster pairs 
or when perfonning null test with randomly re-positioned pairs (see Section ITsl for details). The histograms shows distributions of 40 bootstrapped resamplings 
of the original stacking. When stacking galaxy sets of real cluster pairs, the FL region (red dashed fine) displays a galaxy overdensity above the CI region (blue 
solid line) with a significance of ~ Scr. Galaxy number counts of FL and CI at other redshift slices also show similar high significance detection. In null test of 
stacking randomly re-positioned pairs, we don't observe any overdensity in Fl (blue dotted line) over CI (red dash dot line). Because the galaxy overdensity in 
cluster fields, the galaxy count in Fl and CI from stacking real cluster pairs are higher than stacking randomly re-positioned pairs. See the electronic edition of 
the journal for a color version of this figure. 

luminosity distribution, we first get such distributions of the 

total galaxy population in the FL region and the CI region, 

and then subtract the CI region count from the FL region 

count. For example, to get the color distribution of pure fila- JCL\g~n = 

ment galaxies at one redshift bin, we first get color distribution 

of the FL region with 



f¥L{g-r)-- 



('j)inFL k 

X; 0.49(/riMpc)2 

(ij) in FL 



(13) 



and color distribution of the C 1 region with 



fc\{g-r) = 



E T.gzi^i,yjAMr)kAg-r)) 

('j)inCl k 

E 0.49(/riMpc)2 

(M)inCl 



and then subtract fciig- r) from /fl (g - r) as in, 
Uiiig-r) = fv\.{g-r)-fc\{g-r). 



(14) 



(15) 



/pfl(§-r) is the color distribution of pure filament galaxies 
after removing foreground/background galaxy contents. In 
Equation [T3] and Equation [T4l the 0.49(/r'Mpc)^ factor is in 
the denominators because one pixel in the x and y dimension 
of g(_x,y,M,.,g-r) is 0.7/r'Mpc x 0.7/r'Mpc. This subtrac- 
tion is also illustrated in Figure |5] To acquire properties of 
cluster galaxies for comparison with filament galaxies, sim- 
ilar subtraction between the CL region and the C2 region is 
appUed; 



fciig-r)-- 



Y,ckmr)k:ig-r)) 

_k 

0.49(/2-iMpc)2 ' 

E Egz(Xi,Yj,(Mr)k,(g-r)) 

(U)inC2 k 

E 0.49(/riMpc)2 

(iJ) in C2 



(16) 



/pel (g-r) = /cL ig-r)- fciig - r) . 

Here, fpdig - r) is the color distribution of pure cluster galax- 
ies. 

2.4. Justification on Using Cl/C2for 
Foreground/Background Subtraction 

When using C1/C2 for foreground/background subtrac- 
tion, a concern arises that C1/C2 not only contain fore- 
ground/background galaxies, but also filaments which are not 
in the FL region. In this section, we show that such scenarios 
do exist but the detection of filament galaxies is still effective. 

Indeed, more often than not, galaxy filaments also appear in 
regions other than the FL fields because each cluster is typi- 
cally connected to more than one neighbor In Figure|6la), we 
show the distribution of pair connections for each cluster in 
our cluster pair catalog. Most clusters have more than 5 pair 
connections, and some of them even have more than 10 pair 
connections. However, the probability of one cluster having 
a filament connection in CI region is relatively low. In Fig- 
ure |6tb), we show the distribution of cluster pair connections 
appeared in CI and FL while placing one cluster pair on the 
y = axis in the FL region. About 30% of the clusters do not 
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Fig. 5. — How we get the "pure" color distributions of filament galaxies and cluster galaxies, (a) shows the color distributions of the FL region and the CI 
region, which are almost indistinguishable. After subtracting the CI region distribution from the FL region distribution, we obtain the color distribution of pure 
filament galaxies, which is shown in (b). Subtraction between the CL region and the C2 region is also applied to get the pure color distribution of cluster galaxies, 
(c) shows the color distributions of the CL region galaxies and the C2 region, and (d) shows the pure color distribution of cluster galaxies after subtraction. These 
plots are based on stacking galaxy sets of clusters in the redshift bin [0.14,0.18]. Unless otherwise noted, the eiTor bars in this figure and the rest of the paper 
are single standard deviation eiTors estimated from bootstrapping the galaxy sets being stacked. Because one cluster on average appears 5.7 times in the stack 
(158,897 cluster pairs versus 55,424 clusters), and the error for clusters tend to be underestimated, we enlarge the error estimation of clusters by a factor of \/5.7 
through out the paper. See the electronic edition of the journal for a color version of this figure. 
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Fig. 6. — (a) Distribution of the number of pair connections each cluster have in the cluster pair catalog. Most clusters are connected to more than 5 other 
clusters, (b) Distiibution of the number of cluster pair connections in the CI region (dashed line) and the FL region (solid line). There is no cluster pair connection 
in CI for 30% of the time, and the average number of cluster pair connections in CI is lower than that in FL. Meanwhile, the FL region always contain at least 
one cluster pair connection since one cluster pair is always aligned along v = at the FL side, (c) Cumulative probability distribution of the difference between 
number of cluster pair connections in FL and CL CI is unlikely to have equal or more cluster pairs than FL (~ 35%). The horizontal and vertical dotted lines 
indicate where number of cluster pair connections in FL and C 1 are equal. 



have cluster pair connection appearing in CI, and the number 
of cluster pair connections in CI tend to be smaller than this 
number in FL. In Figure |6lc), we show the cumulative proba- 
bility distribution of the difference between number of cluster 
pair connections in FL and CI. The probability of CI hav- 
ing more or or equal number of cluster pair connections than 
the FL region is less than 40%. On average, the FL region 
contains L09 1 ± 0.007 more cluster pair connections than the 
CI region. Therefore, the FL region always tend to contain 
more filament galaxies than the CI region, and using CI for 
foreground/background subtraction would still leave proper 
galaxy counts similar to stacking pure inter-cluster galaxies. 

Similarly to CI, the C2 region would also contain 
cluster-pair field galaxies. Using C2 for cluster fore- 
ground/background subtraction will then remove cluster-pair 
field galaxy count in addition to foreground and background 
galaxy count. However, because of the very high galaxy over- 
density in clusters, the absence of such a population have neg- 
ligible influence (^ 1%) and won't change the main conclu- 
sions of this paper 

2.5. Null Test with Random Cluster Pair Re-Positioning 

To test the robustness of our algorithm, we randomly trans- 
late the angular coordinates of the cluster pairs and rerun our 
algorithm. We expect to get a stack of galaxies which are 
completely flat with no overdensity either at the "cluster cen- 
ter" or in the "filament region". 

2.5. L Random Cluster Pair Re-Positioning 

For one cluster pair in the 160,954 pre-masked cluster pair 
catalog, we first generate a random point inside the angular 
region 110° < a < 260° and 0° < 5 < 60°. We then trans- 
late the angular coordinates of the two clusters in this cluster 
pair so that the left cluster is laid on the generated random 
point. After randomly re-position every cluster pair in the 
160,954 pre-masked cluster pair catalog, we mask out these 
pairs which fall into regions with dust extinction in r-band 
larger than 0.4, and make a new pair catalog which will be 
referred as the randomly re-positioned pair catalog. This new 
pair catalog and the original cluster pair catalog overlaps sig- 
nificantly in sky coverage, and have similar physical separa- 
tion distributions, as well as similar cluster pair count in each 
redshift bin. 

2.5.2. Galaxy Set Profile of Randomly Re-positioned Pairs 



We run the algorithm described in Section 12.21 and Sec- 
tion 12.31 with the randomly re-positioned pair catalog, and 
get the S(x,y) contour map as shown in Figure |3lb), and the 
distributions of galaxy number count per [/!~'Mpc]^ in the 
FL region and the CI region as shown in Figure |4] We do 
not observe an overdensity either in the filament region or in 
the cluster region when stacking galaxy sets of randomly re- 
positioned pairs. We therefore state that the treatment of sky 
coverage geometry in our algorithm is proper, and the fila- 
mentary as well as cluster overdensity observed in Section l273] 
are real attributes of cluster pair fields. 

3. RESULTS 

3.1. Color Distributions 

This section studies the g-r color evolution of filaments at 
0.10 < z < 0.42. For SDSS filters, despite the 4000A break 
transition from g-band to r-band around z ^ 0.38, practi- 
cally, g-r works almost equally well for separating blue late- 
type galaxies from red early-type galaxies at 0.3 < z < 0.42. 



iHao et all ( 120101) have a detailed discussion on using g-r for 
detecting the red sequence in cluster galaxies out to z = 0.42. 

With Equations [T3l - [T6l we compute the color distributions 
of galaxy filaments, /pfl(g-r), as well as the color distribu- 
tions of galaxy clusters, f^{g-r), at different redshifts. Re- 
sults are shown in Figure |7] The red sequence population is 
observed in both clusters and filaments. Comparing to galaxy 
clusters, galaxy filaments have a bimodal color distribution 
and a larger blue galaxy population. 

More interestingly, galaxy clusters and filaments both show 
evidence of redshift evolution, with a blue cloud galaxy pop- 
ulation increasing steadily from z = 0.1 to z = 0.4. In galaxy 
cluster s, such evolution is obser ve d as th e Butcher-Oemler 
effect ( iButcher & OemlerilI978allbl 11984 . although the ex- 
act nature of it is much debated. Many embra ce this ef- 
fect as testir nony o f hierarchical clustering ( Kauff mannI 1 995 : 
Bauahetal. 1996: E lUngson et alJl200U IMcGee et aljl2009l: 



Li et al. 2009), but ther e are also strong voices that remain 
crit ical (D ressier 198?; ' Small et al.lll998l lAndreon & Ettoril 
ll999i:rAndreon et al.ll2004 120061) . saving that either the error 
of measurements are underestimated, or this effect is no more 
than selection bias. Nevertheless, our data show a very strong 
redshift evolution of the cluster blue fraction and suggest that 
a similar effect also exists in filaments. At redshifts z < 0.2, 
the blue galaxy cloud is only a small bump comparing to the 
red sequence population. Moving to higher redshifts, the blue 
cloud becomes more prominent, taking up the more than half 
of the population at z > 0.28. 

3.1.1. Blue Galaxy Fraction Evolution 

Figure|7]shows that filaments and clusters both contain blue 
and red galaxies. This figure suggests redshift evolution in 
the filament galaxy population. This section quantifies this 
evolution by measuring the blue galaxy fraction in filaments 
and clusters as a function of redshift. 

To measure the blue galaxy fraction, we determine a 
red/blue galaxy cut to separate the blue/red galaxy popula- 
tion. We first take the g-r color of the red galaxy peak in 
filaments and clusters as their red sequence ridgeline, and 
then, i nspired by measurements o f the Butcher-Oemler ef - 
fect in IButcher & OemleJ (11984b : iMargoniner eFaLl (12001b : 
iGoto et alj ( l2003h . we take the g-r value 0.2 mag bluer than 
the red sequence ridgeline as the blue/red galaxy cut. In Fig- 
ure |7] we have this blue/red galaxy cut marked as a vertical 
solid line in both the cluster and the filament color distribu- 
tions. We compute the blue galaxy fraction in filaments R^a 
and clusters R^^\ per redshift bin by evaluating. 



R. 



pfl^ 



R 



pel ■ 



/'''/pfl(g-Od(g-r) 
!U{g-r)d{g-r) ' 

/Vpci(g-r)d(g-r) 
!U{g-r)&{g-r) ■ 



(17) 



Here s^ is the g-r color of the separation between blue and 
red galaxies in this redshift bin. The integrations in the above 
equations are performed with a five-point Newton-Cotes inte- 
gration formula on the binned color distributions. 

In Figure |8] we show this blue galaxy fraction estimation in 
filaments and clusters at different redshifts. We also fit a lin- 
ear evolution of the blue galaxy fraction with redshift (dashed 
Unes) using the MPFITEXY routine (Williams et al. 2010|) 
(dependent on the MPFIT package from iMarkwarda (i2009i) ). 
Our measurement suggest that the blue galaxy fraction in- 
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Fig. 7. — g-r distributions of clusters (a) and filaments (b) in 8 redshift bins. Filaments are bimodal with their color distribution and contain a larger blue 
galaxy bump than clusters. Also, the filament color distribution shows evidence of redshift evolution in the form of a blue gala xy pop ulation that becomes 
increasingly important at higher redshift. The vertical solid lines in these plots mark the blue /red g alaxy color cut used in Section l3.1.ll Note that the galaxy 
population sampled here is somewhat incomplete at 0.38 < z < 0.42 (See discussion in Section l2.2.6t . See the electronic edition of the journal for a color version 
of this figure. 
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Fig . 8 . — Redshift evolution of blue galaxy fraction in clusters (a) and filaments (b). These two figures show box plots of the blue galaxy fraction evaluation in 
clusters and filaments of 40 bootstrapped stacks at each redshift bin. The bottom and top of the boxes represent the lower and upper quartiles, and the whiskers 
represent the minimum and maximum values of the data, or 1 .5 times the quartiles. Data points outside this range are plotted as open circles. We also fit the blue 
galaxy fraction hnearly to the redshift (dashed lines). The blue galaxy population in both clusters and filaments increases at higher redshift, b ut the c hange in 
filaments is less dramatic. Also note that the observed galaxy population is somewhat incomplete at 0.38 < ; < 0.42 (see discussion in Section l2.2.6t . See the 
electronic edition of the journal for a color version of this figure. 
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creases with redshift with a linear fitting slope of 1.18 ±0.05. 
For filaments, the linear fit of the blue galaxy fraction 
as a function of redshift slightly prefers a positive slope, 
0.66 ±0.41, with a significantly higher intersect, which re- 
flects the overall bluer filament populati on. W e briefly discuss 
the implications of this result in Section|4T| 



3.2. Luminosity Function 

In this section, we study the luminosity functions of the fil- 
ament and the cluster galaxy populations. Because our data 
does not allow us to constrain the 3-dimensional distribution 
of cluster and in particular of filament galaxies, we make the 
following assumptions when computing the number density 
of galaxies: 



1 . A^-body simulations predict that the density profiles fall 
off quickly bet ween l/i~^Mpcand2 /i~'Mpc from the in- 
tercluster axis (IColberg et al.ll2005l) . Thus, a single fila- 
ment will typically not occupy the whole FL region. We 
assume that one filament occupies 1/3 of the filament 
region defined in Section |231 and extends ±l/2"'Mpc 
along the line of sight. 



tions: 



$fl(M,) : 



(ij) in FL / 

E 0.49(/i-iMpc)2 ' 

(ij) in FL 

j:cL(M,,(g-r),) 



$cl(M,) = 



0.49(/2-iMpc)2 ' 
E E8z(Xi,Yj,MrA8-r)i) 



$ci(M,): 



(/.»inCl / 



E 0.49(/riMpc)2 

(U)inCl 

E j:8z(Xi,Yj,Mn(8-r)i) 

('j)inC2 / 

E 0.49(/i-iMpc)2 

(iJ) in C2 

I ..^^ $fl(M,)-$ci(M,) 

/>pfl(Mr) = , 

' X 27^; X 2/i-iMpc 



(18) 



^CliMr) = 



Kcl(Mr) ■■ 



3 '^ 20% 
$Cl(M,)-$c2(M,) 



80% 



l^x2/riMpc 



In Figure|9] we plot the luminosity distributions of filaments 
and clusters. We also fit the luminosity d istribution of fila- 
ments and clusters with a lSchechteil (119761) function: 



2. Since most galaxy clusters are at the scale of 1 /i ' Mpc, 
we assume that one cluster lies within the aperture of 
^x-+y- < 0.7/1"' Mpc and also extends ±1/1"' Mpc 
along the line of sight. Note that this does not imply 
a strong triaxiality of the cluster but rather limits the 
choice of cluster galaxies to those found in a cylinder 
fully contained within the cluster. 



3. We assume that the purity of the gmBCG cluster cata- 
log is 80%, whi ch is the lower li mit estimated for the 
gmBCG catalog dHao et al.ll2010l) . 



4. We also assume that on average 20% of the cluster 
pairs in our sample are connected by filaments. This 
again is a lower lim it, estimated from the results of 
IColberg et"aD(l2005h . 



0(M,)dM, = O.41n(lO)0*(lO-°'*<'^'-'^'>)"+'e'"""'"""''dM,. 

(19) 
Because our data can only put weak constraints on the faint 
end of the lu minosity function, we chose to fix a = -1 .2, the 
value fou nd bv|Blanton et al.l (1200 1|) and close to the measure- 
ments of iHovle et alj (12005b and iPopesso et alj (|^05). The 
very bright end of the luminosity distribution tends to depart 
from a Schechter function (Yang et al. 2008). We therefore 
restrict our fits to galaxies withMr-51ogjg/2 >-22.3. We list 
the results of our fits in Table [T] The normalizations, $* and 
r-band characteristic magnitude, M*, of filaments and clusters 
are denoted by 0^ , M*f and c/)*, M*^., respectively. 

In two individual redshift bins at 0.10 < z < 0.14 and 
0.14 < z < 0.18 our measurement of Mj* has large errors. Vi- 
sual inspection of these fits suggests that the fixed faint end 
slope may not be a good representation of the data in these 
bins. Especially in the low redshift bins the data, however, is 
too noisy to leave the faint end slope as a free parameter We 
choose to ignore these two bins in our subsequent discussion 
and analysis. 



The assumptions listed here are clearly only rough approx- 
imations adopted for convenience. We therefore stress that 
they only impact the normalization of the luminosity function 
and not the determination of the characteristic magnitudes as a 
function of environment. Consequently, we caution the reader 
against comparing the measured 3-d galaxy number density 
of filaments and clusters in this paper to measurements ob- 
tained in other ways without taking these approximations into 
account. 

Under the above assumptions, the luminosity distribu- 
tions of filaments and clusters in one redshift bin, denoted 
by (/)pfl(M,) and 0pci(Mr), respectively, are computed from 
gzix,y,Mr,g-r) and cl^{Mr,g-r) with the following equa- 



3.3. Cluster Richness Dependence 

This section studies the dependence of filament overden- 
sity on cluster richness. When stacking galaxy sets of clus- 
ters in the cluster pair catalog, in addition to two redshift bins 
0. 10 < z < 0.26 and 0.26 < z < 0.42, we further bin the clus- 
ter pairs into two richness bins: one bin with both clusters in 
the pair having richness estimation between 8 and 10 and one 
bin with both clusters in the pair having richness larger than 
10. Here, by richness, we mean the GM_SCALED_NGALS 
or GM_NGALS_WEIGHTED richness estimation of the gm- 
BCG cluster ca talog, whichever is recommended for a specific 
galaxy cluster (iHao et al.ll2010l) . 

We evaluate the filament galaxy overdensity of different 



TABLE 1 
Fitted Luminosity Function Parameters 



z 


M, cut/mag 


<^f*/(/r'Mpc-') 


M*f-51og,Q/i/mag 


<P*/(lr'Mpc-^') 


M*,.-51og,gWmag 


[0.10.0.14] 


-24.2 < Mr < 


-18.0 


0.37 ±0.24 


-23.49 ±5.40 


5.18±0.13 


-2 1.544 ±0.033 


[0.14,0.18] 


-24.2 < Mr < 


-18.6 


0.32±0.18 


-24.21 ±5.16 


5.36±0.11 


-21.422 ±0.018 


[0.18,0.22] 


-24.2 < Mr < 


-19.2 


0.78±0.15 


-20.87 ±0.25 


5.21±0.15 


-21.369±0.018 


[0.22, 0.26] 


-24.2 <Mr< 


-19.6 


0.73 ±0.17 


-20.80 ±0.24 


5.01 ±0.07 


-21.317±0.015 


[0.26, 0.30] 


-24.2 <Mr< 


-20.0 


0.65 ±0.12 


-20.76 ±0.20 


4.86 ±0.06 


-2!.229±0.014 


[0.30, 0.34] 


-24.2 < Mr < 


-20.4 


0.41 ±0.11 


-21.09 ±0.39 


4.66 ±0.05 


-21. 193 ±0.011 


[0.34, 0.38] 


-24.2 <Mr< 


-20.6 


0.47 ±0.14 


-20.66 ±0.31 


4.67 ±0.06 


-21.137±0.012 


[0.38, 0.42] 


-24.2 <Mr< 


-21.0 


0.45 ±0.21 


-20.71 ±0.38 


4.78 ±0.07 


-21. 103 ±0.016 
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This table lists the fitted luminosity function parameters as discussed in Section[3J2] The first column gives the redshift range of each redshift bin, and the second 
column gives the magnitude cut at this redshift, i.e., the luminosity function fitting range. The eiTors of the parameters are evaluated through fitting 40 resamplings 
of the galaxy sets being stacked at each redshift bin. Error estimation of clusters is enlarged by a factor of \/5.7 since clusters are on-average oversampled by 
this factor in the stack. 
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Fig. 9. — Lumino sity D istributions of filaments (1st and 3rd rows) and clusters (2nd and 4th rows) and fits to the Schechter function (dashed lines) using 
the MPFIT package |Markwardt 2009). The fitting parameters are listed in Table [T] We do not fit for the M,— 51og[Q/7 < -22.3 bright end of the luminosity 
distributions since a LRG population component is possibly present, and it luminosity distribution tend to deviate from the Schechter function JYang et aU2008l) . 
See the electronic edition of the journal for a color version of this figure. 



bins with a quantity Spg, which is computed with equations: 

(/,;■) in FL k.l 



Sc 



E 0.49(/riMpc)2 

(i.j) in FL 
(!,;■) in CI k,l 



(20) 



E 0.49(/i-'Mpc)2 

{i.j) in CI 
Spfl = SpL-Sci- 

This quantity is the 2-d projected galaxy number count per 



[/z"'Mpc]2 of our detected filament galaxy overdensity. How- 
ever, this quantity should not be misinterpreted as the typical 
surface number density of filament galaxies. Not every cluster 
pair being stacked has a filament connection between them, 
and single filaments typically occupy only a small subset of 
the FL region. Consequently, the real projected galaxy count 
of filaments per unit area can be significantly larger than the 
values given by Epfj. In Table[2l we list the Epg between clus- 
ter pairs at the above richness and redshift binning. The Spfl 
between richer cluster pairs tends to be larger than the Spji be- 
tween poorer cluster pairs. Assuming the geometrical distri- 
bution of filaments does not change with cluster richness, this 
result suggests that richer cluster pairs tend to be connected 
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Epfl [(/!-' Mpc)- 



TABLE 2 
I VERSUS RICHNESS AT TWO REDSHIFT RANGES 



Richness range 0. 10 < ; < 0.26 0.26 < z < 0.42 



8<Richness < 10 
Richness > 10 



0.053 ±0.030 
0.075 ±0.016 



0.012 ±0.020 
0.061 ±0.008 



by higher density filaments. 

3.4. Spatial Distribution 

Galaxy filaments can be warped a nd do not always 
strictly lie along th e intercluster axis CPi mbblet et al.l 120041 : 
IColberg et al.ll2005h . Therefore, significant filamentary over- 
densities still can be seen off the intercluster axis. In this sec- 
tion, we rotate the FL region and the CI region away from 
the intercluster axis as described in Table |3] and compare the 
surface number density of galaxies Spo at two redshift ranges 
[0.10,0.26] and [0.26,0.42]. Our resuh is shown in Figure[l0] 
As we move away from the intercluster axis, the observed fil- 
ament population gradually falls off. At 9 > 45°, the filament 
galaxy overdensity is almost consistent with zero. 

4. SUMMARY AND DISCUSSION 

In this paper we developed an algorithm to detect and 
study the properties of galaxies in filaments statistically. Our 
method makes use of the fact that close pairs of galaxy 
clusters have a large probability of being connected by fila- 
ments. This is a robust prediction of numerical simulations of 
structure formation (e.g., Colberg et al. 2005) and has been 
confirmed observati onally in galaxy redshift surveys (e.g., 
iPimbblet et al.ll2004b . This allows us to study filament galax- 
ies statistically without the need to spectroscopically identify 
them. Our algorithm is applicable to all photometric surveys, 
in which large numbers of galaxy clusters can be found. As 
a result we can study the properties of representative filament 
galaxies out to much higher redshifts than is currently possi- 
ble with wide-field spectroscopic surveys. 

We applied our new method to the gmBCG galaxy cluster 
catalog and the photometric SDSS DR8 BOSS galaxy cata- 
log. Although the projected filament galaxy densities exceed 
that of the surrounding field by only ^ 2%, we can detect 
these relative overdensities at high significance (often ^ 5a) 
in redshift bins of thickness Az = 0.04 at 0. 10 < z < 0.42. The 
limit of this redshift range is governed by the gmBCG clus- 
ter finding algorithm and the small survey volume in the local 
Universe at the low redshift end, and by the flux limit of SDSS 
at the high redshift end. Future deep wide-field surveys like 
the Dark Energy SurvejQ will allow us to push this method to 
higher redshifts. 

As an application of our method, we study the color and lu- 
minosity distributions of galaxy filaments, as well as the de- 
pendence of the filament richness on the richness of the clus- 
ters to which they are connected. Finally, we put limits on the 
warping angles of large-scale structure filaments. 

As verification additional to that already presented in Sec- 
tion |2] we also apply our al gorithm to the Millennium Simu- 
lation dSpringel et al.ll2005b . We present this in the appendix 
and only mention here that such tests support the efficacy 
of our methods and also confirm previous suggestions that 
the Millennium Sim ulation is overpopulated with red galax- 
ies JGuo et alj|201 Ih . 

http : //www . darkener gey survey .org/ 



4.1. Color Evolution 

We find that filament galaxies have a clear bimodal color 
distribution. At redshift 0.10 < z < 0.14, the filament g-r 
color distribution is qualitatively comparable to the u-r color 
distribution of brig ht wall galaxies at z < 0.107 presented in 
iHovle etal] ( l2012h : both being bimodal, and the blue cloud 
galaxy taking up ~ 1/3 of the whole population. 

We find clear indications for a Butcher-Oemler effect in 
galaxy clusters at very high significance. Fitting a linear func- 
tion to the blue galaxy fraction we find that it increases with 
redshift with a slope of 1.18 ±0.05 over the redshift range 
studied. This is comparable to the results of iMargoniner et al] 
( 120011) . who studied the Butcher Oemler effect in 295 Abell 
clusters. Although we chose a slightly larger cluster aperture 
(0.7/r'Mpc instead of 0.49/i~'Mpc) and brighter magnitude 
limit (-24 . 2 < M,. < -20 .4 compared to M*-l KM,- <M* + l, 
which is approximately -22.2 < M, < -19.2 according to our 
measurements of M* ), our value for the slope agrees with their 
result of 1 .24 ± 0.07. Care needs to be exercised when com- 
paring our measurement to some other works studying the 
B utcher-Oemler effect measurements. A later measurement 
by [Hansen eTal] (120091) showed that from z = 0.28 to z = 0.2, 
there is a '^ 5% decrease in the blue galaxy fraction. How- 
ever, our linear fit sugge sts a ~ 9.5% c hange over this redshift 
interval. In their paper, [Hansen et al.l adopted a very differ- 
ent cluster aperture, which depends on cluster richness, and is 
most of the time smaller than our cluster aperture (0 .4/2~' Mpc 
to 0.6/i~' Mpc compared to our choice of 0.7/2"' Mpc). Since 
measurements of the Butcher Oemler effect are likely to be 
dependent on di stance from cluster centers (Ellingso n et alj 
l200UlGoto erariii2003 ; De Pro pris et al. 2004), the difference 
between the measurements of [Hansen et al, and ours is not 
particularly surprising, but a more specified study on this ef- 
fect in clusters may wish to refine the choice of cluster aper- 
ture here. 

We marginally detect an evolution in the blue fraction of 
filament galaxies with a linear slope of 0.66 ±0.41. This 
indicates that the population of blue galaxies in filaments is 
increasing at higher redshift, but this redshift evolution fea- 
ture is weaker than in clusters. The implications of such a 
observations are complicated, though. Should galaxy fila- 
ments be considered as ensemble of galaxy groups and galax- 
ies in less dense environments, our result seems to be sup- 
portive of a Butch e r-Oem le r effect i n galaxy groups, as dis- 
cussed in iLi et alj ( 120121) . iLi et al.l measured the red/blue 
galaxy fraction in galaxy groups at different mass ranges and 
found that galaxy groups of all different masses show a trend 
of increasing blue galaxy fraction at higher redshifts. Their 
measurements of the blue galaxy fraction in galaxy groups at 
their lowest total steller mass bin log(M*,g,p/MQ) < 1 1 .2 are 
-- 48% at z = 0.25 and '- 52% at z = 0.35, somewhat lower 
than but consistent with our interpolated values ^ 57% at 
z = 0.25 and ^ 63% at z = 0.35. Or, if we consider galaxy 
filaments as genuine large scale structures themselves, some 
processes that are thought to cause the Butche r-Oemler ef- 
fect in clusters, such as galaxy galaxy merg ers (iLaverv et alj 
119921: lKauffmannlll995l:lBaugh et alj |199d). should also be 



investigated for filaments. Finally, the observed redshift evo- 
lution of the blue galaxy fraction in filaments is possibly a 
"by-product" of changes in filament galaxy contents or mor- 
phologies at different redshifts. Our results on the evolution of 
the blue galaxy fraction in filaments are currently limited by 
the signal-to-noise ratio per redshift bin. Future deep surveys 



15 



TABLE 3 
Definitions of the FL and CI boxes used in Section 



Number 






FL 






Cl 


1 

2 
3 
4 


2.1/r 
2.1/r 
2.1/r 
2.1/r 


'Mpc<.t < 10.5/1" 
'Mpc<.t< 10.5/!" 
'Mpc<;c< 10.5/1" 
'Mpc<.t< 10.5/1" 


'Mpc, 0° < |e = arctan()'/,v)| < 15° 
'Mpc, 0° < 9 = aicta.n(y/x)\ < 30° 
Mpc, 15° <\e = arctanCy/A:)| < 45° 
Mpc, 30° < 9 = ai-ctanCv/;!:) < 60° 


-10.5/r 
-10.5/r 

-10.5/r 
-io.5r 


'Mpc < a: < -2.1/1" 
'Mpc < a: < -2.1/1" 
'Mpc<A:<-2.1/r 
'Mpc < a: < -2.1/1" 


'Mpc, 0° < |e = arctan(>'/A:)| < 15° 
'Mpc, 0° < e = arctan(y/A:) < 30° 
'Mpc, 15° <\e = arctanCv/jc)! < 45° 
'Mpc, 30° < 6» = arctanCy/^:) < 60° 



0.6 


' 1 ■ ■ ■ 







- XI. 0°^larctan(y/x)l^ 15° 
n2. 15°^larctan(y/x)l^ 30° 
A3. 30°^larctan(y/x)l^ 45° 

. O 4. 45°^larctan(y/x)l^ 60° 
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Fig. 10. — The 2-d projected galaxy number density of filaments, Spfl, with FL and CI boxes placed at different extension angles. The exact definition of FL 
and CI are listed in Table|3] Filament galaxies still can be seen at increasing angles off the intercluster pair axis, but eventually die out when the extension angle 
from this axis is larger than 45 ° . See the electronic edition of the journal for a color version of this figure. 



will give us a longer lever arm in redshift to study this effect 
in more detail. We also expect that improvements in galaxy 
cluster selection will lead to an enhanced catalog purity at all 
redshifts. This will improve the SNR at all redshifts. 

4.2. Luminosity Function 

We measured the characteristic r-band magnitude of clus- 
ter and fila ment galaxies in redshift bins by fitting a standard 
ISchechteil luminosity function wi th a fixed faint end sl ope of 
a = -1.2. For cluster galaxies, iPopesso et al.l (l2005h . who 
studied the RASS-SDSS clusters to z = 0.26, is a suitable com- 
parison to our work. Their measurements of the cluster r-band 
characteristic magnitude at z < 0.26 are -21.35 ±0.19 mag or 
-21.40 ±0.20 mag, depending on the choice of background. 
This is in good agreement with our measurements. An in- 
crease of M* of cluster galaxies with redshift is seen at high 
statistical significance. This result, however, may be affected 
by a potential redshift dependence of the purity of the gmBCG 
catalog, as well as a redshift dependent richness threshold in 
the cluster finder A detailed study of these factors is beyond 
the scope of the present work and the redshift evolution of M* 
should be treated with caution. 

In redshift bins where we obtain reliable measurements, the 
r-band M* of filament galaxies is significantly lower than that 



TABLE 4 
SDSS R-BAND CHARACTERISTIC MAGNITUDE MEASUREMENTS 

M* Measurements 



Filament M* (our measurement) 

Blantonetal. (2001) M.* 

Montero-Dorta & Prada (2009) °'M.* 



-20.85 ±0.11 
-20.83 ± 0.03 
-20.71 ±0.04^ 



" Note that this magnitude is measured for the SDSS r-band blue shifted by 
0.1, i.e., "'r. jBlanton et alJ200"33 IH). and are not directly comparable to the 
rest in the table. At i = 0. 1 , the K-coiTection to such bands ai'e —2.5 log m 1.1 = 
-0.1035, but the K-correction for our system at z = 0.1 ranges from —0.06 
to 0.37. Here, we make a simple assumption that there is a unifonn offset 
between the galaxy absolute magnitude measured in the "'r band and the 
regular r band. Since the the K-con'ection we use at z = 0.1 scatters around 
0.08, which is the K-correction for galaxies of z = 0. 1 and g-r = 0.8 using 
the Chilinsarian et al. (201(|) approximation, we assume the offset between 
the tw o systems to be 0. 1 8. Th e r-band characteristic magnitude measured in 
Monte ro-Dorta & Pradal J2009I) is now closer to our filament measurement. 



of cluster galaxies. We find no significant redshift evolution 



i-.f 



ofM, 

When we stack the luminosity distributions of filament 
galaxies over 0. 18 < z < 0.42, we obtain a characteristic lumi- 
nosity of filament galaxies of M*^ = -20.85 ±0.11 mag. This 
value is compatible with two r-band M* measurements (listed 
in Tablelll) of SDSS galaxies, in which galaxies were not sep- 



16 



arated by their envi r onmen t. This re-enforces the result of 
lAragon-Calvo et all ( 120101) that most matter and thus most 
galaxies live in filaments and not in clusters. Consequently, 
the filament environment is expected to dominate globally av- 
eraged measurements. Our value is also close t o the r-band 
M* measurement of wall galaxies in lHoyle et al.ll2005) , who 
found M* = -20.62 ± 0.08 mag. The filament gal axies in this 
work fall into the category of wall galaxies in iHovle et all 
(l2005h (galaxies in denser environments than extremely un- 
derdense cosmic voids). Not surprisingly, filament galaxies 
are significantly brighter than void galaxi es, whose r-band M * 
was measured to be -19.74 ±0. 1 1 mag by lHoyle et al.l (120051) . 

4.3. Filament Richness and Geometry 

We find that richer clusters are connected to richer fila- 
ments. This i s consistent with d escription of structure forma- 
tion given bv lBond et al.l ( 119961) . In this picture, superclusters 
with pronounced filaments between their massive constituent 
clusters are the result of rare initial density peaks. In turn, 
less massive clusters were "disadvantaged" from the begin- 
ning by being more isolated and having fewer and less dense 
filaments connected to them. Because we select up to 5 neigh- 
bors for any cluster to form cluster pairs with, independent of 
the cluster richness, both factors - the density and the number 
of filaments - affects the measured surface density of galaxies. 
Consistent with this expectation, we find that the correlation 
of filament richness with cluster richness exists at two redshift 
ranges studied here. 

Finally, we also studied the distribution of matter in what 
we call the filament region between the clusters. Slightly 
more than half of all filaments are curved and some connect 
the clusters at which they terminate with offsets from the in- 
tercluster axis. These reasons make it impossible to disentan- 
gle the density profile of filaments from the distribution of fil- 
ament geometries in our stacks. Nevertheless, the galaxy den- 
sity in the filament region should decrease as one moves away 
from the intercluster axis. We verified this by varying the po- 
sition angle of the filament region when rotating it around the 
cluster center We find that filaments are unlikely to be warped 
by more than 45° from the intercluster axis. 

4.4. Outlook 



This stacking approach to studying filaments is not only 
applicable with optical galaxy data, but can also be applied 
to weak lensing measurements, Sunyaev-Zeldovich effect ob- 
servations, and X-ray emission from the warm-hot intercluster 
medium, provided appropriately large survey data exist. We 
look forward to its broader application in the era of combined 
data from deep wide-field optical, SZ, and X-ray surveys. 
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TABLE 5 
Fitted Luminosity Function Parameters in Simulation and SDSS 



4>*[{lr^Mpc ^] (Simulation) 
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0.27±0.13 
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Mi cut 
-25.3 < Mi < -19.5 
-25.3 <M, <-20.1 
-25.3 <Mi< -20.5 
-25.3 < M; < -20.9 
-25.3 <M, <-21.3 
-25.3 <Mi <-21.7 
-25.3 <Mi <-21.9 
-25.3 < Mi < -22.3 



M* 



51ogio/! [mag] (Simulation) 0J[(/)'Mpc"^] (SDSS) M*2-51ogio/i [mag] (SDSS) 



-20.39 ±0.23 
-20.60 ± 0.47 
-20.49 ±0.28 
-20.43 ±0.31 
-20.00 ±0.14 
-19.48±0.31 
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-20.88 ± 0.44 
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[0.18,0.22] 
[0.22, 0.26] 
[0.26, 0.30] 
[0.30, 0.34] 
[0.34, 0.38] 
[0.38, 0.42] 



of the Schechter function fitting in simulation and SDSS perspectively, and M*j and M*, are the characteristic magnitude of the fitting in simulation and SDSS . 
The errors of the parameters are evaluated thi'ough fitting 40 resamphngs of the galaxy sets being stacked at each redshift bin. 
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APPENDIX 

DETECTION OF FILAMENTS IN SIMULATION DATA 

As further test of our method as well as for comparison with simulation data, we run our algorithm on the Millennium sim- 
ulation (Springel et al. 2005), a large volume, high resolution N-body simulation under ACDM cosmology. This simulation 
includes information of dark matter structures, i.e., Friends-of-Friends (FOF) groups, and dark matter substructures. It also 
tracks the merging history of dark matter structures and substructures for galaxy formation simul ation. For our purpose , we 
use the all sky ligh tcone realization of this simulation constr ucted with r nethod s fully described in iHenriques et al.l (120121) and 
iBlaizot et alj ( l2005h . under the galaxy formation simulation of lGuo et alj ( 1201 Ih . 

Due to differences between these simulation data and SDSS observational data, we make a few adjustments of our methods 
before running our algorithm. First, instead of using gmBCG clusters to build the cluster pair catalog, we search for pairs of 
central BCGs in cluster-sized FOF groups with M2()o (the mass within the radius where the enclosed average density is 200 times 
the critical density of the simulation) larger than 1O''*M0. We make this adjustment because running the gmBCG cluster finder 
on these simulations does not yield satisfying purity and completeness, possibly due to the lack of clear color bimodality in the 
simulated galaxy population. Secondly, the simulation galaxy catalog has a apparent magnitude limit in i-band at / = 2 1.0 mag. 
Therefore , when making the brightness cut to ensure fair comparison of the galaxy population at different redshifts (as described 
in Section l2.2.6l ). the luminosity limit is adjusted to -25.3 < M, < -21 .9. This luminosity limit is above the / = 21 .0 mag apparent 
magnitude limit at z = 0.38. For a fair comparison, we also rerun our codes with SDSS data under the new luminosity cut. When 
comparing results from simulation data and SDSS data in this section, all the results presented are under this new luminosity cut. 

In Figure (TT] we sho w the color distri butions of filament galaxies in the Millennium simulations compared to those in the 
SDSS data. As noted by iGuo et al.l (1201 Ih . the simulation seems to contain too many red sequence galaxies in filaments, and the 
fraction of the blue galaxies does not match that seen in SDSS. 

In addition, filament galaxies in the simulation appear to be fainter than observed in SDSS. In Figure[T2] we show the luminosity 
distributions of filament galaxies in simulation and SDSS, and fit them to a Schechter function with the faint end slope fixed to be 
a = -1 .2, as having been discussed in Section [3?2l The fit parameters of the distributions are listed in Table|5] Filament galaxies 
in simulation appear to have a fainter characteristic magnitude in i-band than in SDSS. Possibly because of the magn itude cut 
-25.3 < Mi < -21.9 being too shallow, the Schechter function fittings yield smaller estimations of (f>* than Section [3^ 
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Fig. 1 1 . — g — r distributions of filaments at different redsfiift bins in simulation (a) and SDSS (b). Comp aring to SDSS data, the simulation is over-populated 
with red galaxies. Note that the galaxy population is incomplete at 0.38 < z < 0.42 (See discussion in Section l2.2.6t and also (b) differs from Figure|7](b) because 
of the different luminosity cut applied. See the electronic edition of the journal for a color version of this figure. 
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Fig. 12. — Luminosity distributions of filament galaxies in simulation (1st and 3rd rows) and SDSS (2nd and 4th rows) and their Schechter function fits (dotted 
hues). The Schechter function fit has a fixed faint end slope, and the rest of the fitting parameters are listed in Table |5] The characteristic magnitude of the fitted 
Schechter function is brighter in simulation than in SDSS, indicating overpopulation of bright galaxies in the simulation. See the electronic edition of the journal 
for a color version of this figure. 



